Quantum Charge Fluctuations in a Superconducting Grain 
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We consider charge quantization in a small superconducting grain that is contacted by a normal- 
metal electrode and is controlled by a capacitively coupled gate. At zero temperature and zero 
conductance G between the grain and the electrode, the charge Q as a function of the gate voltage 
V g changes in steps. The step height is e if A < E c , where A and E c are, respectively, the 
superconducting gap and the charging energy of the grain. Quantum charge fluctuations at finite 
conductance remove the discontinuity in the dependence of Q on V g and lead to a finite step width 
oc G 2 A. The resulting shape of the Coulomb blockade staircase is of a novel type. The grain 
charge is a continuous function of V g while the differential capacitance, dQ/dV g , has discontinuities 
at certain values of the gate voltage. We determine analytically the shape of the Coulomb blockade 
staircase also at non-zero temperatures. 

PACS numbers: 73.23.Hk, 74.50.+r, 72.10.Fk, 74.20.Fg 



I. INTRODUCTION 

In the conditions of Coulomb blockade, the charge of 
a grain is well-defined and discrete. It can be varied by 
means of the external parameter, gate voltage V g . Peri- 
odically in V g , the ground state of the system approaches 
a point of degeneracy, in which two consecutive allowed 
values of charge yield the same energy of the ground 
state. In the limit of no tunneling between the grain 
and particle reservoirs (leads), the degeneracy between 
the ground states is indeed reached at the correspond- 
ing special values of V g . However, electron tunneling be- 
tween the grain and leads may remove the degeneracy. 
The nature of the resulting many-body ground state is 
sensitive to the spectrum of excitations in the grain and 
in the leads. Except for the case of ultra-small grains?- 
the level spacing between the single-particle excitations 5 
in metallic grains is usually negligible. Under this condi- 
tion, there are the following three known types of evolu- 
tion of the ground state with the variation of the external 
parameter V g . 

If the grain and leads are normal metals (N-I-N junc- 
tion) the system at the charge degeneracy point can be 
mapped on the multi-channel Kondo problem^ In this 
mapping, two subsequent values of charge and the gate 
voltage play the role of the pseudo-spin and magnetic 
field, respectively. (Normal grain charging was exper- 
imentally studied in Ref. 0-) Similar mapping is also 
possible^ for a system consisting of a normal lead and 
a superconducting grain (S-I-N junction) in the case of 
relatively small charging energy E c < A (A is the super- 
conducting gap in the spectrum of excitations); tunneling 
then lifts the degeneracy between two states with charges 
differing by 2e. The corresponding 2e steps in the grain 
charge were observed experimentally?^ but there was no 
detailed study of the charge dependence on V g near the 
degeneracy point. 

In the opposite case of a fully superconducting system 
(S-I-S junction) with E c < A, the degeneracy removal is 



equivalent to the formation of an avoided crossing in a 
two-level system. Charge degeneracy for such junctions 
was studied experimentally^ 

The third and the last out of the studied classes of 
junctions is represented by a system consisting of a nor- 
mal grain and a superconducting lead (N-I-S junction), 
in which case the degeneracy is not removed by tunnel- 
ing. Indeed, the ground state is degenerate with respect 
to changing the number of electrons in the grain by one, 
whereas the particle reservoir may supply only pairs of 
electrons £ 

We demonstrate here the existence of a new class of be- 
havior of the ground state. It occurs in an S-I-N junction 
having normal lead made of a metal and characterized by 
charging energy E c > A. An adequate model for such 
a junction has large number of channels. We show that 
in the limit of infinite number of channels, the depen- 
dence of charge Q on V g has no discontinuities, but the 
corresponding differential capacitance, C^m — dQ/dV g , 
remains singular, exhibiting a jump at some value of V g . 
We find the full dependence Q(V g ) at a finite (but small, 
G <C 47re 2 /ft) conductance. The smearing of the steps 
in Q(V g ) due to quantum fluctuations may be observed 
in experiments of the type performed in Refs.0 and|^ at 
higher values of the junction conductance, or by using the 
sensitive charge measuring techniques of Ref. y- We also 
find the evolution of the Q(V g ) dependence with temper- 
ature, and find conditions at which thermal fluctuations 
do not mask the quantum effects. 

It may seem that within the framework of the con- 
stant interaction model employed here and in the limit 
of small mean level spacing in the grain the shape of the 
Coulomb blockade staircase in the S-I-N case should not 
differ from that for the N-I-S system studied in Ref. 0. 
Indeed, the charging energy of the system, see Eq. @ 
below, can be reexpressed in a similar form in terms of 
the number of electrons in the lead. Upon this procedure 
the Hamiltonian of the S-I-N system is formally identi- 
cal to that of an N-I-S system. The physical difference 
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between these cases is that, despite the small mean level 
spacing in the grain, for the S-I-N system it is possible 
in the experiments to reach the low temperature regime, 
T < T* = A/ln(V87rA/(5), at which no thermal quasi- 
particles are present in the superconductor i For the N-I- 
S system with a macroscopic superconducting lead such 
regime is beyond experimental reach because the mean 
level spacing in the lead, Si, is many orders of magnitude 
smaller than that in the grain. 

The paper is organized as follows. We present a sim- 
plified derivation of our main results in Section [H] In 
Section [illl we present a rigorous analysis of the problem, 
justify our main approximations and evaluate corrections 
to them. In Section llVI we present a derivation of our fi- 
nite temperature results. In Section[V]we summarize and 
discuss our main results. 



II. QUALITATIVE CONSIDERATIONS AND 
MAIN RESULTS 



In the absence of tunneling between a superconduct- 
ing grain and a normal lead, the system is described by 
Hamiltonian 



H = H C + H s + Hi 



(1) 



where H c , H s , and Hi describe respectively the charging 
energy of the grain, including the dependence on the gate 
voltage, the BCS state in the grain, and the lead, 



H C = E C {N-Mf, H s =J2^k' 



(2) 



Here c pa and c pa are the creation and annihilation oper- 
ators for electrons in the lead, 7^ and -fko- are the cor- 
responding operators for the Bogolubov quasiparticles in 
the superconducting grain; indices p and fc denote or- 
bital states in the lead and grain, respectively, and the 
spin indices take values a = ±. We assume that the elec- 
tron spectrum £ p has a constant density of states near 
the Fermi level (which is a reasonable assumption for a 
metallic lead), while the Bogolubov quasiparticles have 
a gap A in their spectrum, e k = yQ. + A 2 , where 
is the electron spectrum (with the constant density of 
states, too) in the absence of superconductivity. The 
electron number operator in the grain is denoted by N. 
The electrostatic energy of the grain is of the order of 
E c = e 2 /(2C) and depends on the gate voltage V g via 
the term AT — C g V g /e in Hamiltonian H c , where C and 
C g are the total and gate capacitances, respectively. 



A. position of the step 

In the limit of vanishingly small single particle mean 
level spacing in the grain, S — > 0, the ground state energy 



of the system is periodic in gate voltage Af with period 
2. It is also a symmetric function of N — 21 with I be- 
ing an integer. Therefore, when describing the shape of 
the Coulomb blockade steps it suffices to study the steps 
between the charge plateaus with 21 and 21 + 1 electrons 
in the grain. The shape and position of the 21 — 1 — > 21 
steps is then readily obtained using the aforementioned 
symmetry properties. 

Without tunneling, the charge of the grain Q = eN 
commutes with the Hamiltonian and thus is a con- 
served quantity. Minimizing the system energy with re- 
spect to the discrete number of electrons N, we find the 
positions of the steps in charge. The transition between 
the "even" and "odd" plateaus with 21 and 21 + 1 elec- 
trons in the grain, respectively, occurs at 



A/" = 21 + 



E c 



2E e 



(3) 



In the case of a normal grain (A = 0), the steps are lo- 
cated at equally spaced half-integer values of Af. At a 
finite A < E c , the "odd" plateaus of the Coulomb stair- 
case become narrower, the positions of the steps being 
shifted"! by ±A/(2E C ). 

Tunneling between the lead and the grain results in a 
correction to the ground state of the system. The ex- 
citation spectra of the superconducting grain and nor- 
mal lead differ from each other, leading to the difference 
in this "vacuum correction" between the even and odd 
states of the system. To evaluate the correction, we in- 
troduce the tunneling Hamiltonian iJ t , 



Hi 



a ^pa 



h.C. 



(4) 



kpa 



The electron annihilation operators in the grain are re- 
lated to the Bogolubov quasiparticle operators by 7^ = 
Ukdka — Wka*_ k _ a , where — k labels the time-reversed 
state of k, and u k ,v k = (l/y/2)(l ± Cfc/ £ fc) 1/2 - The tun- 
nel matrix elements t kp are related to the conductance of 
the junction, 



G= 2nH 9 ' 9 = 87T 2 J2\hp\ 2 S(C k )S(Q. (5) 

kp 



We are interested in the corrections to the ground state 
energy introduced by H t near a degeneracy point. Such 
"vacuum corrections" appear in the second order of the 
perturbation in H t and result from the formation of vir- 
tual electron-hole pairs across the junction. In the even 
state, the prevailing (because of the charging energy) 
type of pairs has a hole in the normal lead, while in the 
odd state holes are predominantly created in the grain. 
In the latter state, the creation of a pair involves taking 
out an electron from the condensate, which costs extra 
energy 2A, see Fig.^ The resulting vacuum corrections, 
AE e and AE a , at the even and odd sides of the step, 
respectively, are different from each other. This can be 
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FIG. 1: Picture of the virtual electron-hole pairs contribut- 
ing to the "vacuum correction" of the ground state energy 
in a superconducting grain near the charge degeneracy point. 
On the even side of the charge degeneracy points (left panel) , 
without coupling, all electrons in the grain are accommodated 
in the condensate. In the presence of the coupling, the virtual 
states consist of a hole created in the lead and an electron in 
the grain. On the odd side (right panel), one quasiparticle 
is inevitably present above the gap and the prevailing virtual 
states consist in creating a hole in the grain and an electron 
in the lead, with energy cost at least 2A. In comparison, the 
virtual states which result from annihilation of a quasiparticle 
already present in the grain and creating an electron in the 
lead do not have such energy cost. However, by construction 
the number of such virtual states is much smaller than the 
number of states with the energy gap 2A. Therefore the for- 
mer do not contribute significantly to the vacuum correction. 



easily seen in the second-order perturbation theory in H t , 



AE - AE e 
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1 1 



(6) 



E+ + 2E C E+ E_ + 2E C 



where E± = e + £ ± A. The corrected step position, A/"*, 
to the first order in g is defined by equation 



2E c (Af* - Af ) = AE D - AE e 



(7) 



In the limit A <C E c , the result of evaluation of Eq. © 
yields 



AT 



Af + t^-et- 

47T 2 E c 



I + In 



4£ c 



> Afn 



(8) 



At [E c - A) < E c , the factor [1 + ln(4E c /A)] in Eq. © 
should be replaced by v / 21n[l/(3 - 2y/2)] ~ 2.49. Equa- 
tion JHJ should be viewed as the first two terms of the 
pcrturbative expansion for the Af* (g) function. 

When the conductance of the contact increases, the 
odd plateaus become shorter, as it is shown by Eq. 1(8}. 
At the same time, the plateaus acquire a finite slope 
which can also be calculated in the second-order pertur- 
bation theory in H t and is of the order of g/ir 2 at integer 
values of Af. The slope must be small for the plateaus to 
be well defined. When E c — A <C E c , odd plateaus are 



completely suppressed at g exceeding 7.9(1 — A/E c ). De- 
pending on the ratio A/E c , this condition can be realized 
while the even plateaus still remain flat. 

In the following, we neglect the small correction to 
the average grain charge Q related to the finite slope 
of the plateaus, and determine the large (of the order 
of 1) variation of Q which gives the shape of the step 
separating the "even" and "odd" plateaus in the Q vs. 
V g dependence. 



B. shape of the step at zero temperature 

The shape of the step is described by the dependence of 
the average grain charge, Q — e(N), on the gate voltage. 
At zero temperature, it can be found^ by differentiating 
the ground state energy of the system, 



(9) 



The shift of the step position evaluated above comes from 
the grain charge fluctuations with a typical energy of the 
order of A. On the other hand, the shape of the step 
is determined by the low-energy excitations. The band 
width for these excitations is controlled by the closeness 
to the (shifted) charge degeneracy point . This separation 
of energy scales allows us to derive an effective low-energy 
Hamiltonian, Hq, in which the renormalization of the 
step position is already accounted for. The new Hamil- 
tonian Hq acts in a narrow energy band of width D -C A, 
and is designed to describe the low-energy physics of the 
system near the charge degeneracy point, 



2E c \Af-Af*\ < D. 



(10) 



The bandwidth is defined with the energy reference in 
Hamiltonian Hq set as the energy of the odd state ob- 
tained from Hamiltonian H + H t in second order pertur- 
bation theory. That is, energy E c (2l + 1 -A") 2 + A + AE a 
is subtracted from the initial Hamiltonian. All the states 
whose number of electrons on the grain differs from 21 
or 21 + 1 have energy higher than 2E C . Thus, they are 
excluded from the low-energy subspace. Moreover, as a 
consequence of I? -C A, the even states cannot accom- 
modate any excitation in the grain (it would cost at least 
the energy 2A), while the odd states accommodate ex- 
actly one excitation. Starting from Eqs. Q and J3J and 
using the Schrieffer- Wolff transformation, we are able to 
derive the low-energy Hamiltonian, see Section ITTT1 Here 
we present only the relevant for the discussion terms of 
that Hamiltonian: 

H = e |0) <0| + £ H \ka) (ka\ + £ £ pC t ffCpCT 



ka 

P fc p i fccr ^ 



pa 



\c pa +t* kp cl a \0){ka\], (11) 



Here eo is the energy of the grain in the "even" state |0); 
it accounts for the renormalization of the step position 
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by virtual electron-hole pairs with energy exceeding D. 
Up to small terms of the order of g\J AD, this energy is 
e = 2E c {Af - TV*), cf. Eq. (g). The energies of allowed 
states for the Hamiltonian Eq. (|llfl are within the band 
of width ~ 2D. States \ka) have an excitation in state 
(fcer) of the grain with energy ek — A « C 2 /(2A). The 
tunnel matrix elements tk P = u^tkp ~ (l/v2)ifep account 
for the coherence factors values at small energies. Note, 
that the Hamiltonian allows only for zero or one 
additional electron in the grain. 

With the change of energy reference, Eq. © defining 
the average grain charge must be replaced by 



(N) = 21 + 1 



1 dE a 
2E C cW ' 



(12) 



where E is the ground state energy of Hamiltonian 1)110. 

In the zeroth order in tk P , the wave function of an 
even state is a direct product of |0) for the grain and 
some state |/) of the Fermi sea in the lead; the ground 
state is |0) <g> |/o). Tunneling terms in Eq. (fTTJ) modify 
the eigenf unctions. In the lowest order of perturbation 
theory, the wave function acquires the form 



I*) = 



A\0)+J2Pk q „c pa \ka) ®|/} 



(13) 



The second term in parentheses in Eq. I|13(l describes am- 
plitude of a quasiparticlc-holc pair state created due to 
the tunneling; this amplitude is small within the per- 
turbation theory In higher orders of the perturbation 
theory, additional electron-hole pairs may be created in 
the normal lead. However, as we show in the next Sec- 
tion, these additional terms are small if the number 7V c h 
of quantum channels in the junction is large. In the sim- 
plest case, the number of channels is of the order of the 
junction area measured in units of the Fermi wavelength 
in the metallic electrodes (see Sec. 11111 for details), and 
the condition N c h ^ 1 is not restrictive. If it is satisfied, 
the wavefunction Eq. 1|13[) is valid beyond the perturba- 
tion theory. Having this in mind, for now we may use it 
as a trial function for an eigenstate originating from the 
state |/) in the absence of tunneling. Then the ground 
state energy is the lowest value of E which solves the 
equation 



E = e 



E 

kpa 



\ikp\ 9(-£ p ) 



E-Q/(2A) 



(14) 



This lowest-energy solution 
E a {Af) = -2E C 5M 




(15) 



is separated at TV < TV* from the continuum of states 
E > allowed by Eq. I|14|) and can be associated with 




-60 -30 

{N-N')/SN 

FIG. 2: Coulomb staircase at different temperatures T = 
(thick line) and T = 10, 20, 50T q . (N) is the average number 
of electrons in the grain, V g is the gate voltage. The plots 
correspond to the following values of parameters: g = 2, 



T* « 



= 0.7, 
0.1 A, T a 



8tyA/S = 10 . With these parameters, 
» 0.001A, and TV* w 21 + 0.8. The cor- 



rection due to the finite slope of the steps at integer values of 
CgV g /e is not accounted for in the plots. 



a bound state. This state is formed by the quasiparti- 
cle (virtually) populating the grain and the correspond- 
ing hole in the lead. Upon approaching the threshold 
value TV = A/"*, the binding energy of this state vanishes, 
E (TV) oc (TV-TV*) 2 . At TV > TV* the hole is not localized 
near the junction any more; the formerly discrete energy 
level merges with the edge of the continuum spectrum. 
The described qualitative change of the spectrum at the 
even-odd transition is identical to the one in a well-known 
problem of single-particle quantum mechanics. Indeed, 
the same transformation occurs with the spectrum of a 
particle attracted to a three-dimensional well upon the 
gradual reduction of the well depth; particle becomes de- 
localizcd at a certain strength of the potential^ and the 
discrete energy level ceases to exist. 

The number of electrons in the grain at zero tempera- 
ture is obtained from Eqs. (|12|l and 1|15|) : 



(N)o = 2l + l-6(JV* - TV)/ 



TV* - TV 
SAf 



m = i 



i 



(16) 



Eq. (|16f) describes the transition between the "even" 
plateau with 21 electrons in the grain and the "odd" 
plateau with 21 + 1 electrons. The dependence of (N)o 
on TV is shown in Fig. [21 At the degeneracy point 
(TV = TV*), the charge is a continuous function of TV, but 
there is a jump in the differential capacitance, C^is = 
C g d(N)o/dAf. We can define the step position by equa- 
tion (N)o = 21 + 7} and characterize the step width by the 
value of W = (d(N) a /dJ\f)~ 1 at the step position. The 
smearing of the steps of the Coulomb staircase occurs in 
the second-order in g, with a typical width Wo = 166 Af. 
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C. shape of the step at finite temperature 

We turn now to the effect of a finite temperature, T, 
on the average grain charge, (N). Equation (|12[1 can still 
be used, except that E should be replaced with the free 
energy F = —TlnZ, where 



Z = Trexp (-/3H ) 



(17) 



is the partition function of the system^ described by 
Hamiltonian Ijlip . 

To compute Z, we must determine the energies of the 
excited states. For that we once again will be using 
Eq. (|13fl as a trial wave function. The corresponding 
cigenencrgy can be represented in the form E = £ / + £ . 
Here £/ is the energy of the "bare" state |/). In the 
limit of a thermodynamically large lead, £ is a solution 
of the already used Eq. (|14(l . where the step function 
0(£) should be replaced by the Fermi distribution func- 
tion. Similarly to the zero temperature case there are 
two classes of solutions at Af < Af* . In the first class 
there is just one discrete solution whose energy at low 
temperatures is £ = Eo, see discussion in the paragraph 
below Eq. (|19|) . The second class is represented by con- 
tinuum spectrum £ = (f /(2A) — £ p . The full energy for 
a state in the second class can be written in the form 
E = £f + C 2 /(2A), where £/ is the energy of the state 
that differs from |/) by the presence of one hole in state 
p. Since the partition function involves summation over 
all states the difference between £ f and £ / drops out and 
we have, 



f 



f,ka 



(18) 



Factorizing the ./^-independent partition function of the 
lead, Zi oa d = J2f e~ l3£f , and integrating the second term 
over k, we get 



Z = Z } 



lead e 



-f3E 



N eS (T)] , N ctf (T) = 



8ttAT 



S 2 



(19) 

where 6 is the one-electron level spacing in the grain. 
A rigorous derivation of the finite-temperature partition 
function is detailed in Section ITVl 

We would like to note that the replacement of #(£) 
by the Fermi function in Eq. (|14f> leads to a small dif- 
ference between the energy of the discrete state and Eq. 
This finite temperature correction is small in T/Eq. For 
the purpose of evaluating the partition function this cor- 
rection can be ignored at all temperatures because at 
T/E Q > l/lniV e ff(T), see Eq. (O, the partition func- 
tion is already dominated by the contribution from the 
states of the continuum, while T/Eq is still small. 

Using Eq. I|12f) . we obtain the average grain charge at 
finite temperature 



(N) = 21 + 1- 



9(Af* - Af) 



1 + N cS {T)eP E °W 



Af* -Af 
6Af 



(20) 




FIG. 3: Temperature dependence of the step position, 
N*{T), and width Wt- The same values of parameters as 
in fig. [5] is used. 



where Eq(N) and f(x) arc defined in Eqs. i|15fl and IKill . 
respectively. 

The temperature dependence of the Coulomb staircase 
is shown in Fig. [3 At T — > Eq. (|2*D)l reproduces the av- 
erage charge, (N)o, given by Eq. (|16f) . At finite temper- 
ature, the odd plateaus become broader. The position 
of the step, defined by equation (N) = 21 + i, results 
from the competition between the energy of the grain in 
an even state, Eq, and the entropy of the large number 
of odd states, given by N c s(T) in Eq. (fH^I . The step is 
shifted from Af = Af* to Af = Af* (T), 



when 



Af*(T)=Af*-^lnN eS (T) 



(21) 



(22) 



that is when the shift exceeds the zero-temperature width 
of the step, Wo = 16oW. Note that the thermal width 
Wt = T/E c of the charge step at T > T q is smaller than 
that at zero temperature, Wt < Wq. We thus expect 
a strong nonmonotonic temperature dependence of the 
step width with the minimum width occurring at T ~ T q . 
The temperature dependence of the step position and 
width are shown in Fig. [3] 

The junction conductance g does not affect the charge 
steps at T 3> T q . At even larger temperature, T ~ T* = 
A/ ln(-\/87rA/<5), the "odd" plateaus were shown to reach 
the same size as the "even" ones4£ 

In the rest of the paper, we aim at giving a rigorous 
derivation of the main results, Eqs. I|15fl . (|16fl . and (|20fl . 
We'll show that they hold for the experimentally relevant 
case of a wide multichannel junction between the grain 
and the lead. In Section IIIII we derive the low-energy 
effective Hamiltonian <|11[1 ; we demonstrate there that 
the corrections in the number of channels to the ground 
state energy, are small. In Section IIVI we derive the 
partition function for our system. 
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III. EFFECTIVE HAMILTONIAN AND THE 
GRAIN CHARGE AT ZERO TEMPERATURE 

In this Section we derive the effective low-energy 
Hamiltonian (jl 1|) which is used to describe the shape of 
the step of the Coulomb staircase between the even state 
of the grain with a charge Q = (2l)e and the odd state 
with a charge Q = (21 + l)e. We show that it gives the 
same physics as the Hamiltonian of the system, H + H t , 
given by Eqs. and 10} , for a wide multichannel contact 
between the grain and the lead. 

The effective Hamiltonian H is acting in a narrow 
bandwidth D, on a limited set of states which was charac- 
terized in Section [n] only low-energy electrons (or holes) 
may be excited in the lead, and the states of the grain 
are only those with zero or one low-energy excitation, 

|6l,£c(AT -A0 + C 2 /2A<£- 

The states outside the bandwidth D are accounted in Ho 
perturbatively. For the present problem, it is enough to 
calculate such contribution in the second order in Ht- 
Therefore, the appropriate Schricffer- Wolff transforma- 
tion of the initial Hamiltonian isii 

H ps M(H + H t )M - MH t (l - M)—H t M. (23) 

H 

Here, M is the projector operator on the unperturbed 
states whose energy lies within the band D. We obtain 

H =H' + V, (24) 

H> = e P e + E ^P ol l lka P Q + ]T t p cl a c pa 

kcr pa 

y = E [tk P Poilc p<7 p e + i* kv pj valka p 

kper 

+ E E W PP' C U C P'° Pi + E W kk'Poll a lk'aPo- 
p^p'(7 i—e.o k^k' a 

Here, P e and P are projection operators on the even 
and odd states with 21 or 21 + 1 electrons in the grain, 
respectively, t kp = u k t kp , and 

eo = 2E c {No -AT) + E \**pt*!± . (25) 

kper 2A 

The prime here means that only excited states in the 
bandwidth D, (i.e., such as \(l/(2A) - £ p \ < D) arc in- 
cluded. Making the summation on k and p, we obtain 



2E C TV - TV* + 



fJ 



V2AD ^ 



2tt 2 2E r 



(26) 



The degeneracy point TV" defined by condition e = 
would slightly differ from TV* . This difference comes from 
the virtual states with energies within the band D; such 



states were accounted for in the evaluation of A^*, but 
do not contribute to eo. This difference is small at small 
(D -^i A) bandwidth. Note, that by construction of the 
Hamiltonian Ho, we need energy eo at JV = Af* to be 

2 

within the band D. This sets a condition D > #xA, 
which does not contradict our initial assumption DcA. 

Perturbation V in Eq. H24(l has terms describing 
electron-hole pair fluctuations across the junction, and 
terms describing scattering off the junction within the 
grain or lead. The latter terms correspond to second- 
order processes in H t generated by the Schricffer- Wolf 
transformation l|23|l . For instance, the term correspond- 
ing to scattering off the lead when the grain is in the even 
state has probability amplitude 



u pp' = -]C* fe p'* 



kp 



e k -A 2E c + e k + A 



(27) 



Here, the first and second term in parentheses correspond 
to virtual excitations with charge (21 + l)e and (21 — 
l)e in the grain, respectively; the double-prime means 
that only virtual states with excitation energy outside 
the bandwidth are included in the sum, it reduces to 
condition e k — A > D for the first term (we also used 
property t kp = t- k _ p valid for tunneling Hamiltonian 
preserving the time-reversal symmetry) . The amplitudes 
w° p , and w kk i have expressions similar to Eq. l(77|) . We 
will put terms proportional to w pp , and w kk > in V aside 
for a while and return to the discussion of their effect at 
the end of the Section. 

Let us calculate the ground state energy of Ho- At 
Af < Af* , the unperturbed ground state on the "even" 
side, |e) = |0) ® |/o), is the direct product of the BCS 
ground state with 21 electrons in the grain, |0), times 
Fermi sea ground state in the lead, |/o). Its bare energy 
is eo- In the presence of V, we determine its renormalized 
energy, E , with Brillouin-Wigner perturbation theoryiii 
It is the solution of the infinite order equation 



Eo = e + (e| V |e) + ^ 

n^0 



(e|V|n) (n\V\e) 



E 



Eo 

e|y|ni>(m|y|na>- 



(28) 



(n \V\e) 



(Eo - e ni ) • • ■ (Eq - e„ 3 ) 



where \n) are eigenstates of H' , with energy e n . Antici- 
pating that, for typical S-I-N systems with multichannel 
contacts, the series in the Eq. (|28|l can be truncated at 
the term of the second order in V . we obtain: 



Eo 



eo 



E 



\tk P \ 8(—£ p ) 



Zp 



kper ^0 2A 



(29) 



Here again, the prime means that only excited states in 
the bandwidth D are included. This equation is identical 
to Eq. (|14J) . As a result of the summation in Eq. I|29|l . 
the bandwidth D disappears from the equation for the 
ground state energy and we arrive at Eq. (|15f) . 
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The above result holds when higher order in V terms 
can be neglected in Eq. ((29(1 . This is the case for a generic 
tunnel junction between a lead and metallic grain. In- 
deed, typically the area of junction S exceeds signifi- 
cantly the square of the Fermi wavelength in a metal. 
The effective number of quantum channels in the junc- 
tion, iV c h ~ Skp 3> 1 provides us with a parameter al- 
lowing the truncation of series Eq. 1(29(1 in the ballistic 
regime. In a realistic setup, electrons are backscattcrcd 
to the junction from the impurities and the boundaries in 
the grain and in the lead. In a typical situation of a small 
junction to a macroscopic lead the the backscattering of 
electrons from the lead to the junction may be neglected. 
We therefore concentrate on the effects of electron returns 
from the grain to the junction. Due to the finite grain 
size such returns are bound to occur. Quantum interfer- 
ence between returning electron trajectories in the grain 
may lead to a reduction of However, we shall see 

that the corresponding contribution to Eq. remains 
small in the parameter = S/gA -C 1. 

We start with the analysis of higher-order perturbation 
theory terms which involve matrix elements t kp only. Let 
us first consider the "even" side of the transition, at Af < 
Af* . The first term in Eq. that we neglect is 



kk'pp' a 



[Eq — ekp] [Eq — e pp '][E — e k i p ] 



(30) 



where e kp = C 2 /(2A) - £ p and e p > p = e + £ p > - £ p , and 
the prime means that |e/c P |, |e/c' p |, |e p ' p | < D. We define 
the correlation function 

K(t, C, C) = tkptlp'tk'p'ti,^ - £ p ) 

kk'pp' 

xsg-frMc-bwe-cv), (3i) 

which simplifies Eq. I|30|l to 



J 



d^'dCdC 



[e«-£ofe-^o][e C '{-£o] : 



(32) 

with = C 2 / (2 A) — £ and = e + £' — £; the prime 
means |e a |, |e^|, \e Ci \ < D. 

The value of J depends on a concrete model of the 
junction. We consider here a thin homogeneous insulat- 
ing layer separating the grain from lead; the appropriate 
barrier potential is U(r) = UqS(z), where z is the dis- 
tance from the interface. The transmission coefficient T 
for such a barrier depends on the angle of incidence for 
the incoming electron characterized by the normal com- 
ponent k z of its momentum k; in the limit of low barrier 
transparency T(k z ) = (k z /mUo) 2 <C 1. The dimension- 
less conductance of the junction is g = TSkp/Air, where 
kp is the Fermi wavevector and S is the area of the junc- 
tion and T = T(fc^). We can introduce the number 
of channels in the junction by expressing the conduc- 
tance g = TiVch in terms of the angle-averaged transmis- 
sion coefficient, T = T/3. This definition of A^h yields 



Neb = iSkp/A-K. In terms of the matrix elements in the 
tunneling Hamiltonian, the above model of the barrier 
corresponds to 



:fep z <5 2 (k|, -P||). 



(33) 



with v = (mkp /2ir 2 ), V, and Vi being the density of 
states, volume of the grain and volume of the lead, respec- 
tively; m is the effective electron mass. Equation <|3l5|) 
accounts for the conservation of the component kn of the 
electron momentum parallel to the barrier. Inserting now 
Eq. 1)33(1 into 1(31(1 . we obtain the correlation function in 
the ballistic regime 



3 5 2 



87T 4 iV ( 



ch 



(34) 



We can now evaluate Eq. 1(32(1 . It yields the dominant 
contribution in \Eq\: 



Jb = ~ 



9 



8ir 2 N ch 



-A In 



D 



Eq £0 



(35) 



In the vicinity of the degeneracy point, at Af* — Af < SAf, 
the energy Eq is much smaller than eo fx g(AD) 1 / 2 , as 

2 

the bandwidth D satisfies the conditions |jjA <C D <C 
A. Therefore, the logarithmic term ln(D/eo) in Eq. 1(35(1 
is approximately constant and of the order of ln(l/<j). 
As the result, for a grain in the "ballistic" regime, the 
contribution 1(35(1 yields a correction to the grain charge 



SQ 



\n(l/g) SAf 

' N ch Af* - Af 



(36) 



which is small outside the region Af* — Af ~ SAf/Nc^; this 
region is much smaller than 5Af. 

The ballistic estimate Eq. 1(35(1 holds if the virtual ex- 
citation in the grain travels a distance shorter than the 
grain size L and the electron mean free path £; for defi- 
nitencss, wc assume L < I. The length of the path of the 
excitation depends on its typical energy given by |25o|- 
Indeed, the velocity of excitation is v ~ vp(/A, and the 
time of travel is limited by l/|i?o|; here vf is the Fermi 
velocity. The limitation on the excitation path length 
sets the condition for the applicability of the ballistic ap- 
proximation, A|2£o| > Vp/L 2 . The condition is violated 
sufficiently close to the charge degeneracy point, 



Af* — Af • 



gvF 
E C L ' 



(37) 



Here we assumed sufficiently large N^, so that the in- 
terval of defined by Eq. 1(37(1 is shorter than SAf. We 
turn now to the estimate of the fourth-order term l(30|) 
deep inside this interval. 

If an excitation bounces off the grain walls many times, 
it is reasonable to expect chaotization of its motion. This 
prompts one to consider ensemble-averaged obscrvablcs, 
rather than their specific values for a given grain. To 
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this end, we need to express J of Eq. I|30|l in terms of the 
correlation functions of the electron states in the grain. 
We start with representing^ the matrix elements tk p in 
terms of the true electron eigenstates (j>t and x p , 



d 2 xd z (j>* k {r)d z x P (r)\ z=0 . 



(38) 



Here, x is the longitudinal with respect to the barrier 
component of coordinate r. Inserting Eq. i|38[l into (JSJ) 
and averaging independently the states in the grain and 
in the lead, we can express the conductance 



9 = 1T^ I d 2 xdVd z d^ c « (r,r') 
8tHz^ J s 



x d z d z ,Kt.~ {r' ,r)\ z=z , =0 , 



(39) 



in terms of the ensemble-averaged Green's function in the 
grain, 

K C (r, r') = --ImGf (r, r') = V ^ fc (r)0*(r')<J(C - CO, 

k 

(40) 

and similarly in the lead K$ = (— l/7r)ImG| f . As the 
conductance is determined by tunneling events taking 
place on the spatial range I/Zcf close to the junction, 
it is enough to take the Green's functions for half-infinite 
spaces with the appropriate boundary condition that 
they vanish at the interface. Averaging each of them in- 
dependently on the disorder in the metals, we recover the 
result g = TN^ of the ballistic regime. At the same time, 
the correlation function i|31[) is strongly affected by the 
disorder^ Expressing Eq. (|31|l in terms of the Green's 
functions in the metals and ensemble-averaging the prod- 
uct of two Green's functions, the authors of Ref. [lj ob- 
tained 



[iVc(r',r) 



128tt 5 i/ J S 2 

+P C -C'(r,r% =z , =0 . (41) 



Here P w (r, r') describes the evolution of the probability 
density to find an electron at a given point r. For dif- 
fusive motion (under the condition £ <C L), the diffuson 
P u (r,r') obeys the diffusion equation, 



(iw-DVl)P u (r,r') = 6 3 (r-r'), 



(42) 



and D = (vf£)/3 is the diffusion constant in the grain. 
At frequencies u> smaller than the Thouless energy D/L 2 , 
the solution of Eq. i|42[l reaches the universal (zero-mode) 
limit where 

2-7T 

P u (r,r')+P_,(r',r)« 7 ^) (43) 

is independent of the coordinates in the grain. In the 
case of a smaller grain, L < t, Eq. I|42|l does not hold, 
but the universal limit for P aj (r,r') is the samei^, and is 



reached at u> < vp/L. With the help of Eq. (|43() . we may 
now evaluate Eq. ll'{2[l to find 



.9 2 m(l/g) 
64tt 3 



2A 



-E 



(44) 



Here <5 = 1/vV is the electron level spacing in the grain. 
Equation yields (in the region Af* — Af <C SAf) a 
correction to the charge of the order 



ln(l/ ff ) / SAf 



Af* -Af 



8 

5A 



( 45 ) 



Comparing Eqs. l|3^|) and (|43|l . we see that they match 
each other at the gate voltage given by Eq. I|37|) . The lat- 
ter defines the crossover between the ballistic and "zero- 
mode" limits, and corresponds to the path length of the 
virtual excitations of the order of L. The correction 
Eq. (|45|l remains small everywhere except a very narrow 
region around the charge degeneracy point, J\f* — M < 
6N/{Nff/ 2 . 

On the odd side of the charge degeneracy point, at N > 
J\f* , the unperturbed ground state is Pooler \ e )i where 
k is the closest to the Fermi level state. In the zeroth 
in l/iVJk limit, the ground state energy is zero. The 
correction induced in Eq by the perturbation V in the 
Hamiltonian H is found from equation 



En 



E 



\hp\ A(£p) 
Eq — eo — £ p 



(46) 



As the tunneling matrix element involves a single state 
k, its solution should exhibit strong mesoscopic fluctua- 
tions. The typical value, however, is easy to find, 



Eq 



1 

'2^2 



gSHl/g). 



(47) 



Therefore the coupling to the lead induces a small shift of 
Af* which scales proportionally to the one-level spacing 
in the grain and is of the order of gS/E c . Beyond this 
small shift, the ground state energy at Af > Af* is not 
modified by perturbation V . 

Now we return briefly to the effect of the terms pro- 
portional to Wkk' and w l pp , in V. Indeed, they do also 



contribute to the Brillouin-Wigner expansion (|28|l and 
give correction to the ground state energy 1)150. In par- 
ticular, at Af < Af* , the fourth-order- in- H t corrections to 
Eq formed with such terms are given by integrals similar 
to Eq. (|32J). They only differ by the energy ranges of in- 
tegration which account for the virtual states outside the 
bandwidth involved in the evaluation of Wkk' and w ppl . 
The resulting correction is of the same order as, and not 
more singular at Af — > Af* than the contributions we 
have already evaluated, see Eqs. l|3l)|) and P^l . 

To summarize this Section, we have demonstrated that 
in the case of a wide junction (A^h ^S> 1) corrections 
to the ground state energy Eq. (|15|l are parametrically 
small. The smearing of the non-analytical A/"-dependence 
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of the grain charge, see Eq. (|16|) vanishes in the limit of 
zero level spacing S in the grain. Identifying |0) and \ka) 
with P e |<E>) and P jt a \&), respectively, where |3>) is the 
BCS ground state in the grain in the grand-canonical 
ensemble, and discarding the terms proportional to Wkk' 
and w pp , in V, we can put Hamiltonian (|24J) in the form 

CD- 

We note finally that in the opposite limit of a point 
contact, 7V c h = 1, instead of the jump in Cdig we would 
find a smooth crossover— of width 6J\f in the region of 
even-odd transition. 



IV. CHARGE STEPS AT FINITE 
TEMPERATURE 

In this Section, starting from the effective Hamiltonian 
(|llfl for a wide, multichannel contact between the super- 
conducting grain and a normal lead, we derive rigorously 
Eq. (|19|) which gives the partition function of the system 
near the charge degeneracy point Af* . 



now defines a large number of excited states and is im- 
practical to solve for each of them. The prime in the sum 
means that only unperturbed eigenstates in the band- 
width D are included in the sum. Here, let us note 
that Eq. (|29|l is obtained from l|49|) for the particular 
set /o = {n pa } determined by the Fermi function at 
zero temperature, n p(J = #(— £ P ), and with the property 
£fa = • At zero temperature, we were only looking 
for the solution with lower energy, corresponding to the 
ground state. 

The partition function (|17|l which is the sum over the 
full set of eigenstates E\ can be expressed in terms the 
exact density of states, v{E): 



Z 



dEv{E)e 



-BE 



(50) 



where v{E) = ^2 X S(E — E\). The density of state is 
related to the exact Green's function Ge = (E — Hq) -1 
at energy E: 



v{E) 



1. 



-ImTrGV, E + = E + i0 H 



(51) 



A. Density of states 

To compute Z, we must determine the excited states. 

Without coupling, the many-body "even" eigenstates 
of the system without excitation in the grain are denoted 
|0, /) = |0) <g> |/). We recall that |0) is the wavefunction 
for the grain in the even state, while |/) is the wave- 
function of the lead corresponding to a particular set 
/ = {ripo-} of electron occupation numbers for each state 
(per). The energy of the state |0, /) is £f + eo, where 



(48) 



pa 



is the excitation energy for the lead in state |/). The 
many-body "odd" states with one excitation in state (ka) 
in the grain are \ka, /) = \ka) ®\f), their energy is £/ + 
C fc 7(2A). 

At finite coupling, even and odd state hybridize. In 
order to find the eigencrgies, we can still use Brillouin- 
Wigncr perturbation theory like we did for determining 
the ground state energy in Section IlIII Starting from an 
unperturbed eigenstate |0, /), we can write the Brillouin- 
Wigner equation for its energy in the presence of the cou- 
pling by replacing {|e),e ,S } with {\Q,f),£ f + e ,E} in 
Eq. (J2SJ. The solutions of this equation are the exact 
eigenenergies. For a wide junction (in the limits 6 — > 
and A^ c h — ► oo), such equation can still be truncated 
up to second order terms in the perturbation V, like in 
Section IlIII However the Brillouin-Wigner equation gen- 
eralizing Eq. J2HJ) 



E = e Q +£ f +Y, 



kqa 



E-£ } +i q -ejw 



(49) 



We evaluate Tr G in the basis of the unperturbed eigen- 
states |0, /) and \ka,f). Thus we can write the density 
of states 



v(E) = Im 



f ka,f 



(52) 



in terms of the diagonal elements of Ge in this basis: 

G E (0J) = (OJKE-Hoy^OJ), (53a) 
G E (ka,f) = (kaJKE-Har^kaJ). (53b) 

In order to evaluate Ge(0, /), we need to introduce the 
additional matrix elements of Ge 

G E (0,f;0,f) = (OJKE-Hoy^OJ'), (54a) 
G E (kaJ;0,f) = (kaJKE-Ho)-^,/'). (54b) 

They solve the closed set of equations 

6 ff , = (E-S f -e o )G E (0,f;0,f) 
i 

- ^ tkqn qa G E (ka, f q<7 ; 0, /'), (55a) 

kqa 

(E-£ f -(t/(2A))G E (kaJ- : 0,f) 
~ E **,(! - n q a)G E (0, .r; 0, /'), (55b) 



Here we defined Sff> = 1 if all the electron occupation 
numbers in the sets / and /' coincides, otherwise 8 ft/ = 
0. Moreover, the set f qa (respectively f qa ) coincides with 
/, except for the occupation number in state (qa) which 



= 
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is set to n qa = (respectively n qa = 1). Inserting l(55bl 
into l|55a|l . and defining the self-energy 



\^kq\ Tlga 



kqa 



E-£ f +£ q -e k /(2A) 



(56) 



we get the equation 



5f f , =(E-£ f -e - £ B (0, f))G E (0, /; 0, /') (57) 



E 

k.p^q.a 



E-£ f + Z q -$/(2A) 



G E (0, [fc/a] 1 " 7 0, /') 



The second term in the r.h.s of this equation gives a neg- 
ligible contribution in the limit A c h — > oo. (This can be 
shown along the same lines as in Section HTH ) As a result, 
we find 

G E (0, /) = G E (0, /; 0, f) = [E- £ f -e - E B (0, /J]" 1 . 

(58) 

The last term in Eq. 1)25(1 defining eo depends on band- 
width D. By subtracting this term from e and adding 
it to S_e(0, S), we write the Green's function i(5*8|) equiv- 
alently 



G £ (0, f) = [E — £f — wo — E B (0, Z)]" 1 , 
where w = 2E c {J\f — A/"*), and 

£ B (0,/) =E B (0, /) + e - wo- 
Inserting Eqs. 1)25(1 and 1(56(1 into 1(60(1 we obtain 



(59) 



(60) 



kqa 



\tkq\ n 



qa 



E-£ f + Z q -(£/(2A) 

\h q \ 2 e(-£ q ) \ 

Z q -(Z/(2A) ■ 



(61) 



At high energies |£ p |, the electron occupation numbers 
asymptotically behave as the zero temperature Fermi 
function. Therefore, the sum in Eq. ((61(1 is convergent 
and we don't need to specify anymore that only the states 
in the low energy subspace are included in it. 

We can determine G E (ka,f) in a similar way. In the 
same limit A c h — ► oo, we obtain 

G E (ka,f) = g E (ka, f) + g E {koJ) 2 

x$>~ fep | 2 (l-n PCT )Gi3(0;/ p<r )- (62) 
p 

The first term in the r.h.s. of Ea. ((62(l coincides with 
the Green's function for the unperturbed state \ka, /) 
with one excitation in the grain, g E (ka, /) = [E — £ f — 
Cl/{2A)}- 1 . The second term comes from its hybridiza- 
tion with the many-body states with no excitation in the 
grain. 



Inserting Eqs. ((59(1 and 1(62(1 into 1(52(1 . we express the 
density of states as the sum of different contributions: 



u(E) 
8v f (E) 



E 



c 2 

2A' 



-Im 



GeJOJ) 1 



5v f {E) 

dS E+ (0,f) 
dE 



, (63) 
(64) 



In Eq. I(63|) . the sum of 8- functions comes from the first 
term in Eq. ((62(1 while the terms 8vf{E) come from 
Eq. I(59|l and the second term in Eq. I(62|) . 

In the next Section, we use Eqs. 1(63(1 and l(64|) to eval- 
uate the partition function. We show that the contribu- 
tion of 8uf(E), Eq. ((64(1 . to the partition function can 
be approximated by that of a single state with energy 
E = Eq + £ f and derive Eq. ((19(1 for the partition func- 
tion. Finally, we show that Eqs. 1(19(1 and ((20(1 accurately 
account for the thermodynamic quantities of the grain in 
the entire temperature range T < T*, and at any sign of 
Af - N* as well. 



B. Partition function 

Inserting Eq. I(63|) into l(50|) . we obtain the partition 
function 



fOO 



03C2/2A) 



+ / dwe- f}u '8i>f(w + £f) 



(65) 



Introducing the reduced partition function of the grain, 
Z = Z/Z lcad and inserting Eqs. (JSHJ) and (JSUJ into jB3fl 
we obtain 



Z = N QB (T) - 



dto- 



-Im 



1 — d£/(w+)/dw 



,CJ+ - UJQ - E/(W+) / rp 

(66) 

Zi ca d and A ff(T) were defined in Eq. fB^I . In the sec- 
ond term of Eq. I(66|) . (■ ■ -)t denotes thermal averaging 
over with the Hamiltonian of the isolated lead and we 
introduced the reduced self-energy 



E / (w) = E w+£/ (0,/). 



(67) 



The integrand in the second term in Eq. 1(66(1 is a 
regular function of the occupation numbers n p<7 . We 
can expand the integrand in series in the set of n per . 
Then we replace them by their thermal average f(£ p ) = 
[l+exp(/3£ p )] _1 (Their fluctuations {[n pa — /(£ p )] 2 )t scale 
as the inverse number of electrons in the lead and can be 
safely ignored in the thermodynamic limit.) Finally, we 
re-sum the scries and obtain 



Aeff(T) 



oo 



7T 



11 



in terms of the thermally averaged self-energy 



E(w) = (£/(w)>, 



(69) 




\tk q \ 2 0(-^ q )\ 



- C,7(2A) £ 9 ~ 



C fc 2 /(2A) / 



In 169|) we can replace the summations over k and p by 
integrals, then integrate over then take the integrals 
over £ p by parts. Finally, we get 



E( W+ ) = 



9 



16ttT 



V2A(-q; + - Q 
cosh 2 (£/2T) 



(70) 



Equations (|68[) and H70|) express the thermodynamic 
quantities of the system in terms of two definite inte- 
grals. They form the main results of this section. Below 
we show that the grain charge may be approximated by 
Eq. P?J|l. 

If the Boltzmann weight is removed from the inte- 
grand in the second term of Eq. I|68[) it integrates to unity. 
Indeed, at g — 0, this term corresponds to the density of 
states of the single state of the grain without quasiparti- 
cle on it. Therefore, by counting argument it must still 
correspond to a single state at finite coupling. It is easy to 
see that only the spectral weight of the Green's function 
residing at frequencies to < — LU{ = — Tin A ff(T) results 
in non-negligible contribution to the partition function 
in comparison with that of the first large term. We can 
therefore restrict the integration over u> in the second 
term of Eq. tou< uii- 

For uj < — cui the integral in the last expression can be 
readily evaluated and yields 



£(«+) 



5^ 



a) 



■ 9 



AT 
2T 



,0w 



(71) 

It is clear from Eq. (|71|) that the integral over uj in 
Eq. (|68|l converges at the lower limit. Its integrand has a 
pole precisely at Eq < when Af < Af* , see Eq. JOU ■ If 
Eq < —uji the pole E lies within the region of integra- 
tion over ui and gives the dominant contribution to the 
integral which equals to exp(— (3Eq). At higher temper- 
atures, the integral in the second term in Eq. i|68|) may 
not be evaluated by taking the residue at the pole. How- 
ever at these temperatures, the second term is negligible 
in comparison with the first one. This is also true at 
Af > Af* when Eq = 0, and thus we obtain Eq. fT9"|t . 



V. CONCLUSION 

We studied charge quantization in a small supercon- 
ducting grain with charging energy E c > A contacted 
by a normal-metal electrode in the limit when the sin- 
gle particle mean level spacing in the grain, S, is small. 
At zero conductance g of the junction between the grain 
and the electrode the steps in the charge vs. gate voltage 



dependence, Q(A~), are sharp and positioned in the di- 
mensionless gate voltage A~ at Af = 2l±[l/2 + A/(2E c )]- 

At a finite conductance and zero temperature, the 
steps are shifted by an amount oc gA/Ec, Eq. (jHJ), mak- 
ing the odd charge plateaus even shorter. The charge 
steps become asymmetric and acquire a finite width 
Wo oc g 2 A/Ec- We find the shape of the Coulomb block- 
ade staircase, Q(A r ), for the experimentally relevant case 
of a wide junction with a large number of tunneling chan- 
nels, N c ^. In the A c h — * oo limit the ground state en- 
ergy of the system can be determined analytically and 
is given by Eq. (|15|l . The resulting grain charge in the 
ground state is given by Eq. I|16[) . Although the charge 
steps are broadened at g ^ and the dependence Q{M) 
becomes continuous, the differential capacitance remains 
singular, displaying discontinuities at certain values of 
the gate voltage. 

At finite temperature T < T* = A/ ln(V87rA/<5), we 
obtain analytic expressions for the partition function in 
terms of two definite integrals, Eqs. I|68f) and (|70fl . The 
resulting grain charge can be approximated by Eq. J3UJ). 
The shape of the charge step is plotted in Fig. @ for 
several characteristic values of the system parameters. 
At temperatures exceeding the characteristic "quantum" 
temperature scale T q , Eq. I|22|l . the steps acquire an ad- 
ditional thermal shift described by Eq. H21JI . Finite tun- 
neling also leads to a non-monotonic temperature depen- 
dence of the steps width with the minimal width achieved 
at T ~ T q . The temperature dependent shift of the step 
and its width are plotted in Fig. |3 

The existing experiment on such systems^ perhaps 
is not sensitive enough to see the quantum broadening 
evaluated in this paper. Indeed one may think that the 
saturation of the ratio between even and odd plateaus 
at low temperature observed in the experiment was due 
to the quantum fluctuations rather than due to an im- 
purity creating a state within the gap, as suggested by 
the authors of Ref. |(J This explanation however calls for 
fairly large junction conductance of the order of g = 5. 
This value exceeds the experimental estimate of 25k£7 
for the resistance of the junction^ and also would have 
lead to a significant slope of the charge plateaus, con- 
trary to the observations & Nevertheless, the evaluated 
broadening for the N-I-S junction is easier to measure 
than the charge step width in a normal-grain device 
(there the width is exponentially small at small g). Re- 
cent experiment^ mapped out a considerable portion of 
a charge step, broadened by quantum fluctuations, in a 
normal device. We expect that similar experiments with 
a hybrid device considered in this paper may resolve the 
structure of entire step. 
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